Robust transmission of non-Gaussian entanglement over optical fibers 
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We show how the entanglement in a wide range of continuous variable non-Gaussian states can 
be preserved against decoherence for long-range quantum communication through an optical fiber. 
We apply protection via decoherence-free subspaces and quantum dynamical decoupling to this 
end. The latter is implemented by inserting phase shifters at regular intervals A inside the fiber, 
where A is roughly the ratio of the speed of light in the fiber to the bath high-frequency cutoff. 
Detailed estimates of relevant parameters are provided using the boson-boson model of system-bath 
interaction for silica fibers, and A is found to be on the order of a millimeter. 



I. INTRODUCTION 

Gaussian entangled states of two subsystems are well 
studied in the quantum communication and informa- 
tion literature. These states are often encountered in 
quantum communication experiments. However non- 
Gaussian entangled states are also important in the con- 
text of quantum communication. For example, such a 
communication system may be constructed in two steps 
1]: (a) encoding based on product states of Gaussian 
states, (b) decoding based on measurement on this con- 
tinuous set of Gaussian states. During decoding, one 
must generate superpositions of input Gaussian states, 
which are essentially non-Gaussian states. Thus quan- 
tum communication systems may require non-Gaussian 
states. Moreover, it is impossible to achieve a quantum 
speed-up using only harmonic oscillators and correspond- 
ing Gaussian operations Q : the dynamics of such a sys- 
tem can be efficiently simulated classically. This, in turn, 
means that to gain a quantum advantage in this scenario, 
one needs to use non-Gaussian operations. It has also 
been shown that the distillation of entanglement from 
two Gaussian entangled states is impossible using only 
local Gaussian operations and classical communication 
[3- Thus, both quantum speed-up and distillation of en- 
tanglement, which have a close relationship with quan- 
tum error correction, require non-Gaussian operations. 
Recently there have been a few proposals to detect en- 
tanglement in such non-Gaussian states 

In this work we focus on the problem of preserving non- 
Gaussian entanglement in noisy quantum communication 
channels. There have been several proposals based on 
quantum purification protocols and quantum repeaters 
to communicate entangled photonic qubits over long dis- 
tances 0, HI III • Alternatively, entanglement between dis- 
tant nodes can be prepared by measurements along a 
chain of intermediate particles 0, [n|. The problem 
we address here is quite different: rather than setting 
up remote entanglement between distant nodes, we con- 
sider the problem of transmitting entangled field states 
for a long distance along an optical fiber. Specifically, 
we propose a method that preserves the multiphoton en- 
tanglement of a class of non-Gaussian states transmitted 



through an optical fiber. Most optical fibers are known 
to have minimum loss windows at wavelengths of the or- 
der of few microns. The optical frequencies generally 
suffer attenuation inside such a fiber. This loss adds to 
other decoherence processes which degrade the fidelity 
of entanglement transmission. We propose a hybrid ap- 
proach to control decoherence of optical-frequency non- 
Gaussian states inside an optical fiber. Specifically, we 
utilize decoherence-free subspace (DFS) [d| and quan- 
tum dynamical bang-bang (BB) decoupling |l3j protec- 
tion to preserve non-Gaussian entanglement. The DFS 
protection is used against differential dephasing of the 
two field modes used to construct non-Gaussian entan- 
gled states. The BB process is applied in order to deal 
with the remaining relevant decoherence sources, in par- 
ticular Raman scattering. This is done by inserting phase 
shifters at regular intervals along the length of the optical 
fiber, similarly to the proposal in [l4( where this spatial 
BB procedure was used to protect single photon polar- 
ization states transmitted through optical fibers. In this 
manner we provide a novel application of the general hy- 
brid DFS-BB strategy proposed in 0] (see also |l6]). 

The structure of the paper is as follows. In Sec.[n] we 
introduce a class of non-Gaussian states and model their 
interaction with an optical fiber. In Sec. IIIII we describe 
in detail a hybrid approach to eliminating their decoher- 
ence during transmission through an optical fiber. In 
Sec. IIVI we provide a numerical estimate of loss of en- 
tanglement through a realistic fiber using a boson-boson 
model of interaction between the field and the fiber. 



II. INTERACTION MECHANISM OF 
NON-GAUSSIAN ENTANGLED STATES WITH 
FIBER 

A. A class of two-mode non-Gaussian states and 
their entanglement 

The simplest examples of non-Gaussian states of the 
electromagnetic field are the single photon states. Other 
examples are states generated by excitations of a Gaus- 
sian state 0, Another method to produce non- 
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Gaussian states uses state-reduction 0, |2jJ, Hi] . A re- 
cent experimental proposal discusses how to generate 
non-Gaussian states by subtracting a photon from each 
mode of a two-mode squeezed vacuum state j^. We 
consider non-Gaussian entangled states of two field 
modes a and b (a and b are bosonic annihilation opera- 
tors), produced by subtracting p{> 0) photons from one 
of the modes (say, a) of a squeezed vacuum state \<j>), 
given by, 
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where £ is the complex squeezing parameter and P is 
the normalization constant. This state is non-Gaussian 
in the sense that its coordinate representation is a non- 
Gaussian function of space-coordinates x and y, as given 
by, 
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(3) 

where H n (x) is the nth order Hermite polynomial. Note 
that for p = 0, the state becomes Gaussian. 

The entanglement in this state can be verified by the 
Peres-Horodecki criterion [2j|. We find that the density 
matrix of the state under transpose of the b mode 
transforms into 



n,m— 



\n, m+p) (m, n + p\ . (4) 
The eigenvalues of the above matrix can be calculated as 



A - _L|/-|2(«+p) 

^nn — p 2 mi 



(n + p) ! 



(5) 



Existence of negative eigenvalues of the matrix cr reflects 
that the state \4>) is entangled. The negativity AT |24j of 
the state \4>) can be written as 
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(6) 



which is the absolute sum of all negative eigenvalues of 
the density matrix a. The non-zero value of TV reflects 
that the state is entangled. Deviation of the value of JV 
from zero is a measure of the degree of entanglement. It 
is thus clear from Fig. 1 that the two modes become more 
entangled with increasing values of |£|. 




FIG. 1: Variation of negativity N with the squeezing param- 
eter \C,\ for p = 1. 



B. Model of interaction with an optical fiber 

We assume that the state \4>) is transmitted through 
an optical fiber of length L and thereby interacts with 
the fiber modes, which leads to decoherence of the state. 
An optical fiber consists of many dielectric molecules, 
each of which contains many electronic energy levels. 
The interaction of the input field states with the fiber 
can be reasonably described by a molecule-field inter- 
action Hamiltonian. Dominant decoherence mechanisms 
are phase damping and energy exchange between the field 
and the molecules. 

The exchange processes can be expressed via the fol- 
lowing Hamiltonian, under the rotating-wave approxima- 
tion: 



(7) 



where Ai describes the field operators and B>i describes 
the molecular operators. Here the molecules in the fiber 
act as the bath, leading to decoherence. 

The exact form of the field operators Ai in the above 
Hamiltonian depends upon the model of the interaction. 
For example, if the field modes are at near-resonance with 
the single-photon transition in the molecules, dipole cou- 
pling between them leads to A; = a, b. Dipole coupling 
dominates higher order coupling (e.g., quadrupole cou- 
pling or magnetic dipole coupling), which would lead to 
multi-photon processes. Thus, at near-resonance, single- 
photon absorption by the molecules leads to decoherence. 
Clearly, far from single-photon resonance, higher order 
processes can dominate. In light of these considerations, 
we assume that the frequencies of the two field modes 
are so chosen that two-photon processes may occur in 
the system, while the cross-section of single-photon pro- 
cesses becomes negligible. Specifically, we choose the fre- 
quencies of the field-modes to be much smaller than the 
energy gaps between the ground and first excited elec- 
tronic states. 

There are two different kinds of two-photon processes 
that may occur in a molecule: annihilation (creation) 
of two photons described by a 2 , 6 2 , and ab (a 1 * 2 , 6^ 2 , 
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and attf), and photon- number conserving processes (de- 
scribed by alb, aw, a^a, and b^b). If we assume the 
molecules are in their electronic ground states (at low 
temperature, as discussed later), the process of absorp- 
tion of two photons of two orthogonal modes is disal- 
lowed due to certain selection rules [25^ . Moreover, as 
both modes propagate in the same direction through the 
fiber, Doppler shift of the photon frequency causes the 
absorption processes of two photons in the same mode 
to be off- resonance [2^]. Thus, we are led to a situ- 
ation where most of the molecular levels are resonant 
with the second kind of two-photon transitions. This 
means that low-energy scattering of the photons, namely 
Raman (described by a^b and atf) and dephasing pro- 
cesses (described by the number operators n a = a^a and 
n b = b^b, which is equivalent to Rayleigh scattering), are 
most likely to occur inside the fiber. In this case, we can 
write the effective interaction Hamiltonian in the rotating 
wave approximation as 

H, = J29i(<fWi + ab^Bj) + + r b n b )Bi h 

i i 

= Hf ] +Hf\ (8) 

where Bi is the annihilation operator that corresponds to 
relevant transitions of ith molecule in the fiber, g,'s are 
the coupling strengths, B^ h is the dephasing operator for 
the ith molecule (the exact form of this will be discussed 
later), and and T l b are the respective dephasing rates 
of the two modes a and b. 



III. A HYBRID APPROACH TO ELIMINATE 
DECOHERENCE 

A. A decoherence-free subspace against differential 
dephasing 

In this section we discuss a hybrid approach to sup- 
press these two-photon processes, i.e., how the effect of 
the interaction Hamiltonian © on the state \ip) can be 
eliminated. We start with the dephasing processes de- 
scribed by the number operators n a = a^a and n b = b'b. 
We rewrite the dephasing part of the Hamiltonian JSJl, 
J2i(Xa n a + ^l n b)B^ h , in terms of two collective opera- 
tors 

A± = n a ± n b (9) 

as A_ + r + A + )B; h , where P ± = (r* ± T\)/2. 

Note that the two-mode number states \n, n + p) and 
\n + p, n) are eigenstates of the operator A_ = n a — n b 
for all integers p > 0. Thus these states form manifolds 
of decoherence-free subspaces with respect to differential 
dephasing for a given p. In any manifold {\n, n + p}} (p 
fixed), any arbitrary superposition of all possible states 
is also a decoherence-free state 01 • I n this way, in the 
present case the state is a DFS under the action of 



collective dephasing A_. Note that this protection does 
not require encoding, in contrast to the usual construc- 
tion of decoherence-free subspaces . 



B. Bang-bang decoupling of the Raman process 

As a second layer of protection of the state against 
other decoherence processes (described by Raman in- 
teractions), we now follow the general hybrid DFS-BB 
method proposed in In the standard BB decou- 

pling methods 0, one uses very short pulses so as to 
cancel the effective interaction Hamiltonian. Thus one 
ends up with only that component of the total Hamilto- 
nian of system and bath which commutes with the BB 
pulses. If the pulses are appropriately chosen, entan- 
glement generation between the system and bath states 
can thus be prevented, and decoherence of the system is 
prevented. However, one has to apply the pulses in in- 
tervals shorter than the timescale of decay of the bath 
correlation. A recent proposal uses a spatial, rather 
than temporal version of this idea to overcome decoher- 
ence of a single-photon polarization state in optical fibers. 
Ref. shows how to replace the short time-dependent 
pulses with phase shifters at regular intervals. Our ap- 
proach here is similar except that we deal with non- 
Gaussian entanglement rather than single-photon states, 
and introduce a hybrid DFS-BB approach. In this way 
entanglement can be transmitted over a long distance in 
quantum communication systems. 

Let us start with the total system-bath Hamiltonian in 
the form 

H = Hq + Hj, 

H = huj a n a + huj b n b + H B , (10) 

where are the frequencies of the two field modes, Hb 
is the free Hamiltonian of the molecular bath, and Hi is 
given by Eq. (JSJ. The frequencies oj a ^ are chosen so 
that they are off-resonant with single-photon transitions, 
but possibly resonant with Raman transitions in the fiber 
molecules. Here and below we omit the zero-point ener- 
gies when writing oscillator Hamiltonians. 

The initial state is |* ) = \4>) \M), where is the 
input non-Gaussian entangled state and \M) is the state 
of the bath (for simplicity we use a pure-state notation 
also for the bath; below we consider the effect of the 
bath's state more carefully). At the time tl = L/v - 
where L is the length of the fiber, and v is the average 
speed of light in the fiber - the total wave function is 
|*(tx)) = U{(t l ,0) |* ) (the f subscript stands for "free 
evolution"). Here the exact normal-ordered propagator 
is U[(tl,0) —: exp[— i j Q [Hi(x) + Ho{x)]dx] : (in units 
where h = 1). Now, let r = A/v, and let N = L/A > 1 
so that we can expand the propagator as U{(tl,0) ~ 
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iH(NA)r . . . e -iH(2A)T e -iH(A)r^ wnere 

H(kA) = - / [iJo^+i//^)]^ 

A J(fc-1)A 



(11) 



is the average Hamiltonian over the fcth segment. I.e., we 
have neglected deviations from average fiber homogene- 
ity, 4 = ((H(kA) - [ifj(fcA) + H (kA)]) 2 ) (we discuss 
such deviations in Sect ion HVI) . The spatial dependence of 
Hq can come from the bath self-Hamiltonian Hb, while 
that of Hj can come from the coefficients gi and T l a b 
[Eq. ©]. Therefore: 



H (kA) 
H~r(kA) 



hw a n a + frubrib + -ffs(fcA) 

i 

+ J2(^n a + nn b )B l ph 



H^\kA) + H^ , {kA) 



(2), 



(12) 



Moreover, in subsection IIII CI below we argue that Hb 
is effectively a molecular oscillator Hamiltonian, so that 

H B {kA) = J^s^sS, where f2 s = ^ J^ 1 ~ )A fl s (x)dx is 
an average frequency in an s-phonon state. 

To eliminate the Raman processes dynamically, we 
propose use of phase shifters, defined by the following 
operator 



n 



(13) 



which generates a relative phase of it between the two 
modes (alternatively we can define II as e I7r " a or as e mUb ). 
When these phase shifters are incorporated inside the 
fiber at intervals A, the Raman interaction part in the 
Hamiltonian © effectively vanishes. This occurs because 
of the following identities. First, it is simple to show, 
using the Baker-Campbell-Hausdorff (BCH) formula j2(J 
that 



-iit>n a _ i<j> f 



a, (14) 



and similarly for n& and b. Therefore 



n^WT = ( e ^WVe~ l ™° /2 )(e~ l ™ b/2 6e" r " b/2 ) 

= -a f 6, (15) 

Ilatfltf = (e iw W2 ae -iW2)( e -iW2 & t e ^W2) 

= -a6 f . (16) 

These identities imply that the Raman interaction is ef- 
fectively time-reversed every 2 A due to the action of the 
phase shifters. To show the utility of this result note first 
that 



(2) 



(17) 



because the field term of Hq and Hj obviously commute 
with n and because averaging commutes with the phase 



shifter operation. Now, if we install thin phase-shifters 
inside the fiber at positions x = 0, A, 2A, from A to 
B, the evolution will be modified to 

U(t l , 0) pb e - m ( NA > ■ ■ ■ ILe- m(2A *> T Tle- m(A ^Il. (18) 

Note that in writing this expression we have neglected 
the variation of H inside the phase-shifter; this will hold 
provided that the phase-shifter width is much smaller 
than the distance over which deviations 8k from average 
fiber homogeneity become significant. Now assume that 
the average Hamiltonians over two successive segments 
are equal: 



Hf\{k + 1)A) = Hf\kA) 
H~ {{k + l)A) = TTo(fcA). 



(0/ 



1,2 



(19) 



The better this approximation, the better our hybrid 
DFS-BB method will perform; we address deviations in 
Section IIVI In this case, to first order in r, and using 
Eqs. 117|) and (|18fl . we have exact cancellation of H\ be- 
tween successive segments: 

e -{H((k+l)A)T Ue -iH(kA)T U 
e -i7?((fc+l)A)T e -iILff(fcA)nT 

= e ^iT(m(kA) + H { I 1) {kA)+H i I 2) {kA)) 

x e -ir(Ih(kA)-H^ ) (fcA) + H< 2) (feA)) 



-2i(H (kA)+Hy> (fcA))r 



(20) 



where the smallness of t is the justification for adding 
the arguments of the two exponentials in obtaining the 
last line. Using Eq. (|19(l again, this yields the overall 
evolution operator 



U(T L ,0) K e- iTL (H a (L) + HfHL)) _ 



(21) 



Thus, to first order in r (or the inter-phase shifter dis- 
tance A), we have eliminated the Raman term Hj X \ and 

are left with the dephasing term Hj 2 ^ = ^.(r_A_ + 
T + A + )Bp h , as H causes no decoherence. Considering 

(2) 

Hj , we note that the A_ component yields only an 
overall phase on the state In fact, the DFS corre- 
sponding to A_ remains invariant under the action of n 
as n obviously commutes with A_ . Therefore, the overall 
evolution operator reduces to 



U(r L ,0) 



-iTUHo+EiT+A+B^) 



(22) 



(where averages can be considered as being evaluated at 
L) and our remaining task is to eliminate the collective 
dephasing term (proportional to A + ). However, it turns 
out that it is impossible to do this with BB while using 
only linear optical elements. Fundamentally, the reason 
for this is that the group generated by {n a — n^, a^b, b^a} 
(phase shifters and beam splitters) is non-compact, which 
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means that it can at most apply a dilation but not a sign 
change, as required for time reversal in the BB protocol. 
However, as we argue next, the collective dephasing term 
nevertheless has no significant effect on the preservation 
of non-Gaussian entanglement. 



C. Suppression of collective dephasing 

In the low temperature limit, most of the molecules in 
each segment of length A reside in their ground electronic 
states. In thermal equilibrium, the molecular vibrational 
states obey Maxwell-Boltzmann statistics. The state of 
the molecular bath can be well described by a density ma- 
trix p B = Y. s Vs\<i>s){^s\, where p s = ex.p(-E s /k B T)/Z 
is the Boltzmann probability of an s-phonon excitation 
with energy E s — Ml s , ks is the Boltzmann constant, 
and Z — exp(—E s /kT) is the partition function, and 
\4> s ) is the s-phonon excited state. This is equivalent to 
the assumption that the molecular Hamiltonian Hb is 
an oscillator Hamiltonian, i.e., Hb = ^2 s Ml s s (s is the 
phonon number operator), as mentioned in the previous 
subsection. Moreover, we assume that the held couples to 
the molecular vibrations. Then [Hb, -BpiJ = and hence 
Si^phl^s) = s \4>s)- Under these assumptions we have 
that [H , A + Bp h ] = 0, and hence the effect of the pulse 
operator (1221) on the initially uncorrelated system-bath 
state can be written as 



PSB [tl 



0) w U(t l ,0)(\^)(iJ;\®pb)U\t l ,0) 



x(|V>M®l0 a >v&l)e* Ti H °e 



iT L H ir L T + sA^ 



(23) 



Using the expression Q for we hnd, after tracing 
out the bath, that the state of the held modes at the 
hnal time tl becomes 



s n,m 



Xe -irz,[ Wtot +2r +s ](n-m)| n) n + p )^ m+p \ ^ 

(24) 

where Wtot = + ^>b- The hdelity of this state can be 
calculated as 



F(t l ) = (VKtlM 



34 E^ 



pi 



Eici 



l(n+p) ( n + -iT L n(u tot +2r+s) 



which becomes unity if 

tl = ^/{wtot + 2r + s} , 



(25) 



(26) 



where I is an integer. In this expression tl is, of course 
a function of the summation variable s, while in reality 
there is only a single tl- However, we observe that there 
are two physical limits where this dependence of tl on 
s disappears. Namely, if Wtot <C T + then for I — s we 
hnd that r becomes s-independent. In the opposite limit 
of very weak collective dephasing we also hnd that t is 
s-independent (there is a physical upper limit on s in the 
Gibbs state ps = ^2, s Ps\ ( t ) s){ ( t ) s\)- Both limits can be re- 
alized by controlling the held modes uj a ^- Alternatively, 
in the low temperature limit the sum over s in Eq. I|25|) 
is dominated by the vibrational ground state, and again 
the hdelity is unity if t = In/ujtot- 

While the fidelity is generally reduced and recovers its 
initial value of unity only under certain conditions, we 
next show that the non-Gaussian entanglement can be 
preserved provided one knows the value of T+. At the 
exit end of the fiber, the state of the two modes is given 
by p(tl,0) [Eq. JUJ]. Taking the partial transpose over 
the mode b in this state, we obtain: 



s n.m 



xe -iT L [u tot +2r + s](n-m)^ m + p )( mj n + p \ t 

(27) 

The eigenvalues of this matrix are given by 



Ann ~ p2^l ..I 



Xnm = ±-^Vn m (T L T + )\C\ n+m+2p 



(n + p)\(m + »)! , , , . 
rurru 



where we have dehned the phononic visibility factor as 



E 



p s e 



-2ixs(n—m) 



(29) 



For a harmonic oscillator bath this factor can be evalu- 
ated analytically, using E s = Ml(s + 1/2): 



}_ e -hn/2k B T 



1 

2Z 



E« 



-[HQ/ k B T+2ix(n-m)]s 



sinh(/iil/2fcsT + ix(n — to)) 



(30) 



which is an oscillatory function of its argument x, with 
period n/(n — to). 

The negativity at the fiber's end becomes: 



1 oo oo 

N = JT2 E E v nm (r L T + )\C\ n+m+2p 



n— m^n 



(n +p)\(m +p)\ 



n\m\ 



(31) 
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Comparing this expression to Eq. JJ]J , it is apparent that 
the negativity is reduced due to the phononic visibility 
factor. The condition for this factor to become unity is: 



tt/T 



+ > 



(32) 



which can be satisfied provided one knows the value of 
r_|_. One way to extract this value is, in fact, to apply our 
bang-bang protocol and to test the fidelity of the state via 
quantum state tomography (see |27j for a more general 
method relating BB to tomography). Also note that in 
the low temperature limit, where the sum over s involves 
only a small number of terms, the phononic visibility 
factor will still be close to unity. On the other hand, it 
is clear that if neither condition is met (r ^ tt/T + and 
high temperature) then the rapidly oscillating terms in 
the phononic visibility factor will destructively interfere 
and cause entanglement loss. 

Note also that in absence of any BB control, the ini- 
tial state \ip) evolves through the interaction Hamiltonian 
JSJ). This leads to different combinations of the photon 
numbers in the two modes due to energy exchange with 
the molecular bath. Thus the initial entanglement is de- 
stroyed at a length scale A c , corresponding to the dissi- 
pation time-scale t c of the molecular bath. 



IV. EFFECT OF FIBER INHOMOGENEITY ON 
DECOHERENCE 

So far we have assumed that the average Hamiltoni- 
ans over two successive segments of length A are equal 
[Eq. 119|) ]. However, due to inhomogeneity inside the 
fiber, the average Hamiltonian differs between segments. 
In this section, we show that this fluctuation of the aver- 
age Hamiltonian leads to dissipation of the field and thus 
sets an upper limit to the value of A. We follow and im- 
prove the method described in Appendix A of |14| . 



A. A Gaussian fluctuations model 

The inhomogeneity in the fiber may arise due to 
nonuniform number density of the molecules or slow 
time-dependence of the fiber properties. In view of this, 
we modify the assumption of homogeneity to read 



(I = a, b) are the fluctuations in the coupling coefficients 
in the fcth segment, and e <C 1 is a proportionality con- 
stant defining the strength of the fluctuations. We as- 
sume that e ~ r. These fluctuations lead to losses inside 
the fiber. However, in an amorphous silica fiber, the 
loss due to Rayleigh scattering (i.e., due to terms con- 
taining n a and n&) is much greater than the loss due to 
other mechanisms (e.g., due to Raman scattering, i.e., 
due to terms containing alb and atf) [2^|. With this in 
mind, we neglect the fluctuations in the Raman terms 
and henceforth consider only the effect of the fluctua- 
tions in (r* Thus, in the interaction picture with 
respect to the energy term H + 5Hb (t) , the interaction 
Hamiltonian corresponding to the fluctuations in the kth 
segment of the fiber can be written as 



m)n b ) 



(35) 



where f{(t) = 5r|(t)S!, h (I = a,b) arc the operator- 
valued fluctuations for the ith molecule, and we have 
replaced the /c-dependence with time-dependence, in the 
joint limit r — ► and total number N of segments large, 
such that tl = Nt is finite. Under the action of this 
Hamiltonian SHj(t), the evolution of the field state can 
be described, in the Born approximation, by the follow- 
ing equation pjj ]: 

p F = -iT± B [8Hi{t),p F {Q)®p B {Vij\ 

- Tr B [ [SH I (t),[5H I {t') 1 p F (t')®p B {0)M36) 
Jo 

where the terms up to the second order of e(~ r) have 
been considered (as e <C 1). Now note that since A«L, 
the number of segments of length A is much larger than 
unity. In this limit, the fluctuations fl(t) can be con- 
sidered as described by Gaussian operators . Then it 
is reasonable to assume that the two-time correlations of 
the form Y{t, t') = ( /«(t) £j f> (f)) ) fl (I, V = a, b) 

between these fluctuations do not depend on a, b due to 
the symmetry properties of Gaussian operators, while the 
mean (£\ //(£)) b vanishes. Here (X)b = T^b(pb(0)X) 
denotes the average of any operator X over the bath. 
Taking this average over Eq. I|36|l leads to the vanishing 
of the term linear in SHi(t), while using the expression 
for \ip) [Eq. J3>] we obtain from the integral term: 



H((k-l)A)=H(kA)+d{H B )k + S(H I ) k , (33) p F = - 

where the operator-valued fluctuations S(H B )k,S(Hi)i c 
in the Hamiltonian are independent. Here<5(i/s)fc is the 
fluctuation in the bath-only Hamiltonian, and S(Hj)k 

6(Hj) k = Y, S 9H^bB t + atfB}) k 

i 

+ eJ2(SKn a + STin b ) k (Bi h ) k (34) 

i 

is the operator-valued correction to the interaction 



P 2 



dt'r(i,t')y^C I1+?, C* (m+p) 



n.m 



(n +p)!(m +p)\ / s2i \ / i 
; — ; in — m) n 5 n + p>(m,m + 3?L 

'777.' 

(37) 



n\m\ 



Hamiltonian in the kth segment, where 5g!? and S(T 



l)k 



This can be solved for the matrix elements of p F , yielding 

(PF )n,n+ P ;m,m+ P (r L ) = exp[-4e 2 (n - to) 2 / dt dt' 

Jo Jo 

r(£j t')}(p F ) n ,n+p;m,m+p(G) , 

(38) 
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where the initial state is 



(PF) (0) 



p2 



+p >*(m+p) 



(n + p)\(m + p)\ 

11 1 m I 



(39) 

Let us assume that the molecular bath is in thermal 
equilibrium at temperature T. If T — > 0, most of the 
molecules reside in the ground electronic states. How- 
ever, molecules can be distributed in all the vibronic 
modes corresponding to the ground states. The degener- 
acy of these vibronic states can be lifted by phonon ab- 
sorption due to molecular collisions. In view of this, we 
treat the molecules as bosons inside a fiber. As implied 
by our discussion in the previous section, the dephas- 
ing operator J2i Bph can be written as (B^B)^, where 
(-B)vib is the annihilation operator corresponding to these 
vibronic states. Thus, the interaction of the field modes 
with the molecules can be described by the so-called in- 
dependent oscillator (10) model [3fJ of boson-boson in- 
teractions. In this model, each of the oscillators of the 
passive molecular bath is linearly coupled to the system 
oscillator. This interaction is governed by two terms: (i) 
The Gaussian fluctuations Y^ifii^) °f the bath opera- 
tors, and (ii) a memory function fj,(t) of the bath opera- 
tors, that vanishes at negative times. If the correlation 
between fluctuations at different times vanishes, the in- 
teraction reduces to a Markovian process. However, in 
general, this correlation is a non-trivial function of time 
and thus corresponds to a non-Markovian process. It can 
be shown that the symmetric autocorrelation of // (t) 
satisfies I3CI: 



J(E//(*)E^(*')+E//(*')E//(<) 



ui coth 



2k B T 
Re[jj,(uj + iO + )]duj , 



cos[u(t-t')] 



(40) 



while correlations of an odd number of factors of i ff (t) 
vanish. Here jl(z) is the Fourier transform of the 
memory function /i(t). Due to the Gaussian prop- 
erty of Ylifi(t)> as discussed before, the above equals 
d2ifi(t)12jfi'(t')}- The memory function is indepen- 
dent of the potential and the properties of the system 
and only depends upon the coupling strengths of the field 
operators with the bath. In the 10 model, the coupling 
strength is an even function of the oscillator frequency 
and is of the form irijU) 2 , where rrij is the mass of the jth 
oscillator and the u)j is the frequency of the jth oscillator 
mode. The spectral distribution of the memory function 
is in the Ohmic class and can be written as 



Re[£(w + i0 + )] « mJ > e- u / Ue E '> 



(41) 



Here uj c is the cut-off frequency of the molecular bath 
modes, and is introduced such that for large frequencies, 



the memory function does not blow up. This is in confor- 
mity with the passivity condition of the molecular bath, 
which also requires that the bath modes must have an 
infinite spectrum and the memory function must not be 
a singular function of u> 30]. Thus using Eqs. I|40l) and 
1)41(1 . we can write the following expression for the rate 
of dissipation at time Tl- 



T(t l ) = [ dt [ dt'T{t,t') 
Jo Jo 

duw 3 e~ u/ ^ coth 



2k B T 



2 . 2 f UJT L\ 



(42) 



where we have normalized the memory function (|41ll in 
units of mass. For low temperature (T — ► 0), when the 
quantum fluctuations dominate the thermal fluctuations, 
one can find the actual loss figure as 



T(t l ) « u. 



2 x 2 {3 + x 2 ) 

(1 + X 2 ) 2 



■, X = UJ C T L 



(43) 



Using the results for the unitary evolution Q24J1 and the 
decay of the density matrix elements given by (|38|l . we 
find the following expression for the density matrix of the 
field modes at time tl'- 



p{TL)dh 



52 E^ E c n+p c* (m+p) J ( n+ p)^ m +p) 



p2 



x e -iT L (uj tot +2r + s)(n-m) e -4e 2 r(T L )(n-rn) 2 



1 71, 71 + p) (to, TO + p\ 



(44) 



The eigenvalues of the matrix obtained by taking trans- 
pose of the b mode in the above density matrix are given 
by 



A, 
A„ 



i|CI 2 ^^ v„. 



± 



i 

p2 



p s e 



— 2iTL^+s(n—m) 



ICI 



n+m+2p 



x j {n+p)\( m +p)\ er ^ T[TL){n _ m? ^ 
V n\m\ 

(45) 

Thus in the low temperature limit, we find the following 
expression for the negativity: 



A/- dISS = i E ier +m+2 ^/ ( " +p)!(m+p)! 



n!m! 



xe 



-4e 2 r(Ti)(n-m) 2 



(46) 



where T(tl) is given by Eq. I|43|) . Clearly the negativity 
decreases as the field propagates through the fiber. Us- 
ing the expression l|43|) for T(tl), we find that for x ^> 1 
(for tl S> 1/^ c , i-e., on a time-scale much larger than 
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the bath correlation time), the negativity does not de- 
crease further and saturates to ffify with T(t l ) replaced 
by ui^. This saturation is due to the fact that 1/lo c is 
the timescale over which information about the system 
state spreads in the bath; for times much longer than this 
the bath is effectively stationary and no more damage to 
entanglement in the system is possible via the bath. We 
next extract a distance scale for the phase shifter sepa- 
rations from these considerations. 



B. Numerical estimates of the inter phase-shifter 
distance 

The result (|43[l due to fluctuation given by (|34|) leads to 
an estimate of the spatial separation A between two suc- 
cessive phase-shifters. Note that for n = m (i.e., for di- 
agonal terms of the density matrix) this fluctuation does 
not lead to any dissipation (or loss of entanglement), as 
is clear from <|38|) . But the larger is \n — m\, the higher is 
this dissipation. In other words, the largest off-diagonal 
contribution to the negativity comes from the terms with 
\n — m\ = 1, which therefore are sufficient to give us the 
desired estimate of A. If we allow an error probability 
S(tl) over the length of the fiber, then a sufficient con- 
dition for the present result to be useful is 

e -A^T{r L ) > i_ 6 (r L ). (47) 
Rewriting the above using (|43|l . we have 

x ^^-H1-8{t l )]-V\ (48) 

where we have used e = r. The limiting value is attain- 
able by fixing A and letting the number of phase shifters 
N become very large (long fiber), with uj c given and fixed. 
This sets an upper bound on the applicable value of A, 
which is essentially the ratio of the speed of light in the 
fiber to the bath high-frequency cutoff. 

In the following, we present a numerical estimate in a 
realistic situation, e.g., for an optical fiber with an amor- 
phous silica core. The inhomogeneity in silica leads to 
fluctuation as described above and thus decoheres the 
input fields. 

The Debye temperature = hw c /kB of crystalline sil- 
ica is 342 K, where u c is the maximum phonon frequency 
(frequency "cut-off") allowed inside the crystal. Thus, 
the lifetime r c of phonons becomes of the order of 1/uj c . 
On the other hand, in amorphous solids the Debye tem- 
perature and the lifetime of phonons are not well defined. 
However, at low temperatures Tc6, there exist certain 
empirical relations between them 0,H2. For example, 
at T = 0.2 K, which corresponds to a phonon frequency 
of 2.62 x 10 10 Hz, the life-time of phonons is of the order of 
1CP 10 s. We consider this frequency of phonons as max- 
imum frequency lu c allowed inside the fiber at T — 0.2K. 



3.0-, 
2.8-'; 
55 2.6- 
2.4 - \ 



H.O 2.0 4.0 6.0 8.0 10.0 

CO X. 
c L 

FIG. 2: Variation of the negativity Adiss with time in units 
of the bath high-frequency cut-off, uj c tl for a fixed squeezing 
parameter |£| = 0.5 and p = 1 accounting for (dashed line) 
and not accounting for (solid line) the bath fluctuations. The 
other parameters are lu c = 2.62 x 10 10 Hz, r = 4.325 x 10~ 12 
s, and error probability S — 5 x 10 -2 . 



It has been shown that long distance distribution of en- 
tangled states of two qubits over a noisy quantum channel 
can be achieved using entanglement purification proto- 
cols 0, H| • These protocols can be improved in terms of 
the requirement of physical resources as well as the er- 
ror threshold, if one uses quantum repeaters [jj- It has 
been shown that for an error probability ^0.01 (which 
is much larger than the error threshold for fault-tolerant 
computation using single qubits) inside the communica- 
tion channel, quantum purification protocols work well, 
when combined with quantum repeaters. Although, as 
explained in the introduction, the setup of our problem 
is quite different from that of quantum repeaters, for the 
sake of concreteness we use the threshold figure from that 
scenario and conservatively consider the case when the 
maximum error probability S(tl) allowed through the 
fiber is 5 x 10~ 2 . 

We consider a multimode fiber of length L — I km. 
The time of propagation of the fields through the fiber is 
T£ = 5.33 x 1CP 6 s, where v = c/n g , n g = 1.6 being the 
effective group index of the field through the fiber. Us- 
ing the parameters discussed above, we find from Eq. I|48|) 
A < 0.8 mm. This means that the field sees the phase 
shifters at a time interval r of 4.325 x 10~ 12 sec which, 
as required for the BB protocol, is much smaller than 
the time-scale for bath dissipation, t c = 10~ 10 sec. We 
show in Fig. 2 how the negativity (|46|l varies inside the 
fiber for the parameters discussed here. We find that the 
negativity becomes constant after a certain length scale 
inside the fiber, as discussed before. This is because, in 
the presence of bang-bang control, the effective contribu- 
tion of bath fluctuations to the negativity vanishes at a 
time-scale when the bath correlation vanishes. 
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V. CONCLUSIONS 

In conclusion, we have discussed in detail, how one can 
preserve entanglement in a class of continuous variable 
non-Gaussian states against decoherence caused by cou- 
pling to a bosonic bath. Specifically, we have considered 
transmission of an entangled state of two bosonic modes 
through an optical fiber and developed a hybrid approach 
combing decoherence-free subspaces and bang-bang con- 
trol to sustain the entanglement. We described the non- 



Markovian interaction with the bosonic bath consisting 
of molecules in the fiber and provided a detailed estimate 
of the relevant parameters to implement our approach in 
a realistic fiber. It turns out that to achieve a loss figure 
of < 5 x 10~ 2 in a 1 km fiber, phase shifters should be 
placed about 1 mm apart. This appears to be a tech- 
nologically feasible requirement. Hence we expect that 
the method proposed here will become a useful tool in 
the effort to transmit non-Gaussian entangled states over 
optical fibers. 
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